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Abstract 

General conditions for the occurrence of mesoscopic phase fluctuations in condensed 
matter are considered. The description of different thermodynamic phases, which 
coexist as a mixture of mesoscopically separated regions, is based on the theory of 
heterophase fluctuations. The spaces of states, typical of the related phases, are char- 
acterized by weighted Hilhert spaces. Several models illustrate the main features of 
heterophase condensed matter. 
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1 Introduction 



In many cases, condensed matter does not consist of a unique thermodynamic phase but 
is separated into mesoscopic regions filled by different phases. A number of examples of 
such matter are described in review [1]. The coexisting phases can be distinguished by their 
order parameters that change from one mesoscopic region to another. The order parameters 
are often associated with some special symmetry [2]. However, this is not compulsory, 
and in many cases different phases cannot be characterized by different kinds of symmetry. 
For instance, crystalline and fiuid phases have different spatial symmetry. Magnetic and 
paramagnetic phases possess different magnetic symmetry, which can be observed by means 
of neutron scattering [3]. Similarly, ferroelectric and paraelectric phases possess different 
symmetry characterized by polarization. If the latter phases coexist, the polarization of 
the mixed system is spatially nonuniform [4,5]. Contrary to this, liquid and gas cannot be 
distinguished by their symmetry, though they may have different densities or other properties 
[6]. Metal and nonmetal are also distinguished not by symmetry but by their conductivity. 
Better and worse conducting mesoscopic regions can coexist in many materials [7-9]. 

Mesoscopic is a key word here, implying the following. Suppose that inside one thermo- 
dynamic phase there appear the nuclei of another phase. Each nucleus, in order to represent 
a coherent germ of a phase, should have a characteristic coherence length Icoh much longer 
than the mean interparticle distance a. As is evident, an object of the size a cannot be 
associated with a thermodynamic phase. At the same time, the coherence length is to be 
essentially shorter than the characteristic length L of the whole system. Thus, the phase 
germs are of the mesoscopic size [1] , that is 

a <^ Icoh L . 

This property is to be understood on the average. Generally, the shape of heterophase germs 
is not necessarily regular, but can be rather complicated. Several germs can also touch each 
other forming very ramified fractal-type objects, similar to fractal animals at criticality [10]. 
The total effective length of such a fractal formation can reach the length L of the whole 
system. Being mesoscopic in size, the heterophase nuclei are randomly distributed in space, 
forming no regular structures. 

Random is another key word. Heterophase germs are random in space, analogously 
to random disordering defects in sohds [11,12], but being different from the latter by the 
mesoscopic size of the phase germs. 

These germs, being mesoscopic in space, can also be mesoscopic in time. This means that 
they may be not static but dynamical, arising and then disappearing, with a characteristic 
lifetime, or coherence time tcoh, such that 

Tloc -C tcoh -C 

where tioc is the time of local equilibrium and Tobs is an observation time [1]. The temporal 
behaviour of these fluctuational dynamic germs is chaotic. When the system is chaotic in 
time and random in space, it is termed turbulent. Since here the spatio-temporal chaoticity 
is connected with the dynamics of heterophase nuclei, this phenomenon may be called het- 
erophase turbulence [13]. The origin of heterophase turbulence can be classified into three 
types, extrinsic, intrinsic, and stochastic [13]. 
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Extrinsic origin of turbulence assumes the existence of sufficiently strong external forces 
acting on the system. Such forces can induce chaoticity in classical and quasiclassical systems 
[14,15]. 

Intrinsic origin presupposes that the system behaves chaotically owing to its own internal 
properties. There are many examples of such chaotic dynamical systems [16,17]. For statis- 
tical systems this means that quasiequilibrium heterophase states are their chaotic attractors 
[1,13,18]. 

Stochastic origin of heterophase turbulence is intermediate between the above two. It 
takes into account that all real systems are never completely isolated but can be only 
quasi-isolated, and that an infinitesimally small external stochastic perturbation can trig- 
ger chaoticity, making the system stochastically unstable [1,18,19]. 

The most interesting cases are those related to the intrinsic and stochastic origins, when 
heterophase states are not forced by external sources but arc self-organized. There exist 
plenty of examples of such systems, reviewed in Ref. 1. Among those heterophase materials 
that have attracted much attention in recent days, it is possible to mention liquid crystals 
[20], high-temperature superconductors [21-24] and colossal magnetoresistant materials [25]. 
Mesoscopic phase fluctuations lead to the occurrence of many unusual features of matter, for 
instance such as colossal magnetoresistence [25] and anomalous diffusion in solids [26]. Their 
appearance may transform a second-order phase transition into a first-order one [1]. They can 
induce different anomalies in the behaviour of thermodynamic and dynamic characteristics of 
matter, for instance: the maxima of specific heat below the critical temperature [27] , specific- 
heat jumps at the nucleation point [1,28], anomalous saggings of the Mossbaucr factor at 
the phase transition point [29,30], anomalous decrease of sound velocity and increase of 
isothermal compressibility at the point of a structural transition [29,30], and other anomalies 

[!]• 

All unusual features of condensed matter with mesoscopic phase fluctuations can be self- 

consistently described in the frame of the statistical theory of heterophase fluctuations [1]. 
In order not to force readers to study the whole review [1], a brief survey of this theory is 
presented in Sections 2 and 3. The main purpose of this presentation is, from one side, to 
give a rather short account of the principal ideas, the theory is based on, and from another 
side, to clearly point at the principal mathematical techniques providing the foundation. An 
attempt of combining brevity with accuracy in outlining the theory is a methodical novelty 
of the paper. Several examples of sections 4 to 6 illustrate applications of the theory. The 
generality of the approach to treating mesoscopic phase fluctuations as well as the generality 
of the phenomenon itself are emphasized in Section 7. 



2 Description of Phase Separation 

When talking about phase separation, one, flrst of all, keeps in mind the separation of 
thermodynamic phases in real space. Different phases can be distinguished by their different 
properties. Those phase characteristics that are the most appropriate for distinguishing 
different phases are called the order parameters. The order parameter of a thermodynamic 
phase is a set rj of quantities typical of the phase and uniquely distinguishing it from other 
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phases. If there are several phases, each of them is characterized by the corresponding order 
parameter rj^, with the index z/ = 1, 2, . . . enumerating the phases. When a physical system, 
occupying a real-space volume V, is phase separated, then this volume can be divided into 
subvolumes Vj, filled by different phases with different order parameters rj^. A family {V^} 
of subvolumes Vj, forms an orthogonal covering of V, such that 

v = Uv., v^nv. = vv.. (1) 

The latter can be represented as the resolution of unity 

^e.(r) = l (rev) 

for the manifold indicator functions 

A family of indicator functions 

e = {ar)k = l,2,...; reV} (3) 

uniquely defines a phase configuration in the real-space volume V. 

In physical systems, the separation of thermodynamic phases is usually realized so that 
these phases are not sharply split apart, but are transformed one into another through an 
interphase layer, where the order parameter continuously changes from the value typical of 
one phase to that characteristic of another phase. How could we then locate a separating 
surface between such phases? First, the location of the separating surface can be fixed in a 
unique manner by employing the Gibbs method of separating surfaces [31]. This can be done 
by invoking some additional conditions. For instance, the equimolecular separating surface is 
uniquely defined by requiring that the surface density of matter be zero. Moreover, a precise 
location of the separating surface is not important. This is because the phase configuration 
in what follows will be varied, with the family of indicator functions (3) being treated as a 
random variable. 

In order to define a method of varying phase configurations, we should consider a manifold 
of all possible families (3) and to rig this manifold with a differential measure. To this 
end, we introduce an orthogonal subcovering {y^i} of the volume Vj,, so that 

rii, 

V, = U V,, , V^, n V,,- = 5^, % V,, . (4) 
1=1 

To each subregion Vj^j, we ascribe a vector a,^j e Yi,i playing the role of a local center 
of coordinates, so that moving SL^i implies a congruent motion of V^^. The shapes of the 
subvolumes Nyi as well as the location of the centers a,^j inside each V^i can be arbitrary. 
For a given subcovering (4), the indicator functions (2) can be decomposed as 

C^^W =I]Ci^i(r-a^i) , (5) 

i=l 
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where 



1, reV^i 
0, r^Y^i. 



The differential measure on the manifold {^} is introduced [1] as the differential functional 
measure 

]im 6 - 1 n^^^ nn ^ > (6) 



in which 

X 



, = 1 / e.(r) dr , V = [ dr. 



The measure (6), with varying Xi, e [0, 1] and aj/j e V, induces a topology on the manifold 
{^}. This results in the topological configuration space 

X ^ {i\ Vi} , (7) 

composed of all possible phase configurations. 

There exists another principal problem connected with phase separation. This is how 
such a macroscopic quantity as a thermodynamic phase is related to the microscopic states 
of a physical system. For generality, we consider this question for quantum systems. 

Suppose that all admissible microscopic quantum states of a physical system form a 
Hilbert space 7i. Let be an orthonormal basis in this space, such that ((^m, Vn) — ^mn- 

Then the space of microstates can be presented as a closed linear envelope 

n = Z{ipn} (8) 

over the basis {(fn}- Any state f eH can be written as a decomposition 

f = Y1 f-^^n , fn = {<Pn: f) ■ 

n 

For any two states /, /' e 7i, their scalar product is given by the equation 

(/) / ) fnfn ■ 

n 

Now assume that for each representative of the basis {(fn} a weight Puifn) is assigned, 
normalized so that 

n 

This defines the weighted basis {fmPuifn)}- The weighted scalar product is introduced as 

(/,/').^EP^(^n)/:/;. (9) 
n 

Then a closed linear envelope over the weighted basis results [1] in the weighted Hilbert space 

n„ = C{ipn,Pu{(Pn)} , (10) 
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with the associated scalar product (9). 
Prom definition (9), one has 

For an operator A^,, given on 0.,^, the trace operator is defined as 

I^W. = ^^^n),. ■ (12) 

n 

Prom equahties (11), it follows that 
Therefore the trace (12) reads 

TvH, A^ = J2P''i^ri){^n,A„ipn) ■ (13) 
n 

Using the above formulas, we may specify the weighted average 

<A,> = Trn^ A, (14) 
of an operator Ai,, with pi, being a statistical operator on Tiu, normalized so that 

Trw, Pz. = 1 . 

The order parameters rju can be presented as the averages 

r}v = <fiu> (15) 

of the corresponding order operators fji,. Each order parameter is assumed to possess the 
properties classifying the related thermodynamic phase. These properties implicitly deter- 
mine the appropriate basis weight ((/?„). The latter is not uniquely defined by Eq. (15), 
but, generally, there can exist a class of such basis weights satisfying the same Eq. (15). Any 
representative of this class is acceptable, provided all of them yield the coinciding averages. 
In practice, it is even not necessary to find an explicit form of Pi^{(pn), but it is sufficient 
to impose the appropriate restrictions on the averages, which would guarantee the typical 
features of the related phases [1]. There are several ways of selecting pure phases, such as 
the Bogolubov method of infinitesimal sources [32], the method of thermodynamic quasiav- 
erages [1,33], the method of order indices [34,35], and others [1]. A very convenient way is by 
imposing additional conditions on density matrices or Green functions [36]. A similar way 
of imposing conditions on distributions or correlation functions is also suitable for classical 
systems. 



6 



3 Heterophase Statistical State 

The total space of microstates for a system with phase separation is the weighted fiber space 

y = ®.H, , (16) 

whose fibering yields the weighted Hilbert spaces (10). 

Let Ai,{^i,) be an operator on Tiu, when a phase configuration (3) is fixed. One may also 
consider the operator density Ai,{^^, r) for which 

A,{i,) = jA,{C.,r)dr. (17) 

Here and in what follows, the integration with respect to r, if not specified, implies the 
integration over the whole system volume V. The operator density on the total space (16) 
is defined as 

^(e,r) = ©,^,(e.,r), (18) 

so that the related operator is 

A{0= J A{i,r)dr = (BuA,{^u) . (19) 

Examples of operator densities are the Hamiltonian density H{^,r) and the number-of- 
particle operator density A'"(^,r). 

Statistical properties of a heterophase physical system, with a fixed phase configuration, 
can be characterized by the quasiequilibrium Gibbs ensemble [1]. A particular phase config- 
uration ^ is connected with the local inverse temperature /?(^, r) and local chemical potential 
r). The latter enter in the quasi- Hamiltonian 

Q(0 = / m, r) [H{^, r) - /.(C, r)N{^, r)] dr . (20) 

This defines the quasiequilibrium thermodynamic potential, which will be called quasipoten- 
tial 

q = -^]nTry J e'^^^ , (21) 

with N being the total number of particles and where the integration with respect to ^ is 
over the configuration space (7) of all possible phase configurations. The quasiequiblibrium 
statistical operator is 

h(n = exp{-g(0} 

P^^' Try J exp{-Q{Om ' ^ ^ 

Expectation values of operators A{^) are given by the averages 

< i > = Tr^ / piOAiO Vi . (23) 



Integrating with respect to phase configurations, we assume that measuring procedures 
effectively average the observable quantities over random phase configurations. If all parts 
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of the heterophase system ate locally equivalent to each other, then the local inverse tem- 
perature and local chemical potential, being averaged over phase configurations, should not 
depend on r, which can be formulated as the condition of heterophase equilibrium 

(3^1 r)V^, ii = j /.(e, r) . (24) 

Accomplishing the integration over the configuration space (7), we keep in mind the 
existence of the thermodynamic limit 

N 

N ^ oo , V ^ oo , — — >■ const . 
Then the thermodynamic quasipotential (21) reduces [1] to 

q ^ - 1 InTry e-^^ , (25) 

with the effective Hamiltonian 

H = ®,H, , = H,{w,) . (26) 

The terms 



Hu = J [H^{wu,r) - iJ,N^{wu,r)] dr 



of the Hamiltonian (26) can be called [1] the phase-replica Hamiltonians. The operator densi- 
ties A^{w^, r) have the same form as ^1/(^1/, r), but with replaced by w^,. The quantities w^, 
are the phase probabilities, playing the role of average geometric probabilities for the corre- 
sponding phases. These probabilities are defined [1] as the minimizers of the thermodynamic 
quasipotential (25), that is, 

q = abs mm q{{wi,}) , (27) 

under the condition 

^w,^l, 0<w,<l. (28) 

V 

For the operator averages (23), we find [1] 

< i > = Tr^j; p i , (29) 
where the effective statistical operator 

exp(-/3g)^ (30) 
Tr3;exp(-/3i/) 

is expressed through the Hamiltonian (26). The operator A in the right-hand side of Eq. 
(29) is 

A = ®^A^ , A^ = A^{w^) . (31) 
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Prom Eqs. (26) and (30) it follows that 

Tin. exp{-(3H^) 

Therefore the average (29) can be written as 

<i> = ^ <A,> . (33) 
1/ 

The set {A} of all operators, corresponding to observable quantities, composes an algebra 
of local observables A. The heterophase statistical state is the family < A>. 

The operators from the algebra of local observables can be divided into single-particle, 
two-particle, and, generally, fc-particle operators, so that the overall structure of an operator 

is 

= E ^ ^1'^ ■ (34) 

The standard Hamiltonians of statistical systems usually contain the single-particle and the 
two-particle terms only. Hence the structure of a replica Hamiltonian is 

= Hi'^ + ^wl Hl^^ . (35) 

Assuming that i?^*^^ do not depend on w,^, we may explicitly derive the conditions of het- 
erophase stabihty. 

Consider, for concreteness, a heterophase mixture of two phases. It is convenient to 
introduce the notation 

w = Wi , W2 = 1 — w . (36) 

Then the phase probability w is a minimizer of the thermodynamic quasipotential q = q{w), 
which implies the equations 

^ ^ — 
dw ' dw"^ 

The first of the latter, taking into account form (25), yields 



1^ = 0, H>0. (37) 



The second of Eqs. (37) has the meaning of the positive definiteness of the heterophase 
susceptibility 

X-^>0. (39) 

Taking account of Eq. (25), we have 

2 



^ N 



dw'^ \ dw I 



(40) 
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Then the stabihty condition (39) reads 



> 



(41) 



Introducing the averages 



1 1 

K = < H^^^ > $ = < H^'^^ > 



(42) 



we find 



1 9^/^ 1 0^ H 

— < 7^ > i^i + - ^2 - (1 - W)^2 , < > = ^1 + ^2 ■ 

N dw N ow^ 

The extremum condition (38) yields an equation for the phase probabihty 

$2 + - 
w = — . 

$1 + ^2 

Prom the inequahties < w < 1 it follows that 

-$i <Ki-K2< *2. 

And the stabihty condition (41) gives 



(43) 



(44) 



$1 + $2 > 



I 

N 



\ dw j ^ 



> . 



(45) 



Inequalities (44) and (45) can be called the conditions of heterophase stability. These are to 
be complimented by the usual conditions of thermodynamic stability, such as the positivity 
of specific heat and of isothermal compressibility. 



4 Theory of Melting and Crystallization 

An important application of the developed approach is to describing the phase transitions 
of melting and crystaUization [37,38]. The mesoscopic fiuctuations inside a crystal are the 
liquid-hke droplets presented by the regions of local disorder of the crystaUine lattice, aris- 
ing below melting temperature. Above the latter, mesoscopic heterophase fiuctuations are 
quasicrystalline clusters inside liquid. 

To consider this phase transition, let us take the replica Hamiltonian in the standard 
form 

f ( h^V^ \ 
H^ = w^ J V'i(r) ( - - l^j ^i/(r) dr + 

+1 wl j ^l{v)^l{v')^r - r')^.(r')^.(r) dv dr' , (46) 
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in which are field operators. The chemical potential /i can be defined from the equation 

N^Ni + N2, N, = - <—^> . (47) 

OfJ, 

Here, we deal with a two-phase mixture, because of which u = 1,2. 
The local density of particles in a u-phase is 

p,(r) = < i^i{r)i^.{r) > . (48) 

The number of particles in the u-phase, according to the definition in Eq. (47) , becomes 

^ j p,{v) dv . (49) 

The average density of the z/-phase is 

P^^=l^j P^i^) . (50) 

And the mean density of particles in the whole system writes 

N , , 

p = — = wipi + W2P2 • (51) 

If the density p is fixed, then the latter equation defines the chemical potential p — p{p, T) 
as a function of density and temperature. 

Let the crystalline phase be numbered with z/ = 1, while the liquid phase, with v — 2. 
By this assumption, the local density (48) for the crystalline phase is periodic, 

pi(r + a) = Pa (r) 7^ const , (52) 

where a is any lattice vector. But the liquid phase has a uniform density 

P2(r) = P2(0) = const . (53) 

Since it is the local density that distinguishes the phases, pui^) can be treated as the order 
parameter. Alternatively, it is possible to define the order parameter 

T]^ = sup 1 . (54) 

r Pu 

For an ideal crystal, supi.pi(r) = pi(a), but for a liquid, P2(r) = const. Therefore 

771 = ^ - 1 > , 772 = . (55) 
Pi 

The analysis [37,38] of the typical behaviour of the phase probabilities w^, results in the 
following figTircs. In Fig. 1, the probability of the crystalline phase is shown as a function 
of temperature 6 in arbitrary units. At zero temperature, the system is in pure crystalline 
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phase. For higher temperatures, the crystaUine-phase probabihty shghtly diminishes along 
the hne AB and then abruptly falls down from the point B to E at the melting temperature 
Gm- Above the latter, the curve EF corresponds to the probability of crystalline clusters 
inside the liquid phase. The branches BC and DE are related to metastable states of an 
overheated crystal and crystaUine clusters inside an overcooled hquid, respectively. 

Figure 2 displays the probability of the liquid phase W2 = I — Wi. At zero tempera- 
ture, W2 = 0. With increasing temperature, there appear a small probability of liquid-like 
disordered regions in the crystalline phase, which rises along the line ab. At the melting 
temperature 0^, the probability of the liquid phase jumps up to the point e, after which it 
follows the hne ef. The branches be and de again are related to metastable states. 

In the similar way, one may treat a liquid-glass transition. A glassy state can be described 
by a disordered lattice. Then the local density Pi(r) is not periodic, but, at the same time, 
it is not constant. The glassy solid state and liquid state can be distinguished by the order 
parameter (54), so that, similarly to Eqs. (55), one has 

771 > , 772 = . (56) 

Possible behaviour of the glassy-phase probability is shown in Fig. 3. Contrary to the crystal- 
liquid transition, the glass-liquid transition can occur as a continuous transition, without 
the jumps of phase probabilities at the glassyfication temperature O^. The glassy state is 
metastable, being not an absolute extremum of a thermodynamic potential, as compared to 
the crystalline state. 

Thermodynamic potentials for heterophase matter are the same as for pure-phase sys- 
tems. We can consider, for instance, the grand potential 



the free energy 

or the Kramers potential 



fl = -pV = n{T,V,ii) , 



F^Q + fxN^ F{T, V, N) , 



All of them can be expressed through the quasipotential (25) by means of the equations for 
the grand potential 

n = NqT , (57) 



free energy 

or for the Kramers potential 



F = A(gT + /i), (58) 



K = -Nq . (59) 



As is seen from the latter relation, the quasipotential (25) is minus the Kramers potential 
per particle. From the above equations, the pressure is 

P = - ^ = -PQT , (60) 
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where p is the mean density (51). The minimum of the quasipotential (25) corresponds 
to that of the grand potential (67) or to the minimum of the free energy (58), but to the 
maximum of the Kramers potential (59). 

The quasipotential (25), due to the form of the effective Hamiltonian (26), can be written 
as a sum 

g = , 5, = - 1 In IV^^ e-^^'' . (61) 

Thence, the grand potential (57) becomes 

Q = ^ , = -Tin IV^, e-^^-' . (62) 

And the free energy (58) also writes as a sum 

F = Y,Fu, F, = Q, + pN,. (63) 

V 

However, one should not forget that the terms in the sums (61) to (63) depend on thermo- 
dynamic parameters as well as on the phase probabilities, e.g., 

^l, = n,{T, V, II, w,) , = F,{T, V, iV, w,) , (64) 

where Wy is a function of T, V. fi or T, V. N, respectively. These terms correspond to phase 
replicas of a heterophase system, but not to pure phases. The latter would be described by 
the potentials 

= n'y{T, K, /x) = n^iT, K, /X, 1) , Fl^ FliT, K, Ny) = Fy{T, K, N,, 1) ■ (65) 

These do not compose the thermodynamic potentials (62) and (63) of a heterophase system, 
that is, 

^^E^'u, F^Y.K- (66) 

The difference 

nir^t = ^-Y.K^ Frnt = F-Y.Fl (67) 

V V 

between the total thermodynamic potentials (62) or (63) and the sums of the related poten- 
tials (65), representing pure phases, can be associated with the thermodynamic potentials of 
interphase layers [1]. Analogous conclusions concern the Kramers potential (59), for which 

K^Y.^^n.e-^'''' ■ (68) 

V 

But this potential is a function K{f3,V,f3/j) of not the most convenient thermodynamic 
variables, because of which it is rarely employed in practice. 

In this way, the whole heterophase system is characterized by a total thermodynamic 
potential, like the grand potential (62) or free energy (63). Such a potential writes as a sum 
of the terms representing phase replicas of heterophase matter, so that these terms should 
not be confused with the potentials of pure phases. A phase-transition point, defined in 
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the theory of heterophase fluctuations, can be rather different from a transition point given 
by the conventional approach deahng solely with pure systems. In the latter case, a phase 
transition between two phases, say the melting-crystallization transition, is defined by the 
equalities 

Q[{T, V, //) = Q',{T, V, f,) , F[{T, V, N) = F^{T, V, N) (69) 

for the potentials of pure phases, which give, e.g., the transition temperature T'(/i) or T'{p), 
respectively. This temperature coincides with the transition temperature Tq in the theory 
of heterophase matter only if wi = 1 and W2 = below Tq, while wi = 0, W2 = 1 above Tq. 
However, generally, ^ 0, 1, and the transition point in heterophase matter is different 
from that in pure-phase systems. Thus, the phase transition temperature Tq in heterophase 
matter is defined as the point where the relation between wi and W2 changes so that 

wi >W2 (T = To - 0) , wi< W2 (T = To + 0) . (70) 

In addition, the existence of mesoscopic heterophase fiuctuations may lead to various pre- 
transitional effects [1], which can be explained only taking these fiuctuations into account. 

5 Turbulent versus Chaotic Crystals 

The majority of crystals are associated with a periodic crystalline structure. There exist 
as well quasicrystals, characterized by a quasiperiodic spatial distribution of particles [39]. 
Crystals and quasicrystals can form stable matter. The structure that is not uniform and 
neither periodic nor quasiperiodic, can be called chaotic. Examples of solids with chaotic 
structure are glasses and amorphous materials. In nature, such chaotic solids are mctastable. 
Ruelle [40] put forward a question asking if solids with random spatial structure could form 
not metastable but stable states. Newell and Pomeau [41] mentioned that such static chaotic 
spatial structures could exist in principle, which was supported by a two-dimensional numer- 
ical modeling [42]. The solids with chaotic spatial structures can be named chaotic crystals. 
Sometimes, they are also termed as turbulent, which, however, is clearly a misnomer. Turbu- 
lence is a spatio-temporal phenomenon characterizing unstable temporal evolution of chaotic 
spatial structures [43]. It is, therefore, incorrect to call a static spatial structure turbulent. 
To dignify a phenomenon by the name of turbulence compulsory implies the involvement of 
chaotic temporal evolution. Therefore stable solids with a static random spatial structure 
[40-42] can be correctly classified only as chaotic crystals. 

The existence of temporal heterophase fiuctuations is what can render solids really 
turbulent. Randomly fiuctuating in time and chaotically distributed in space mesoscopic 
heterophase fiuctuations realize heterophase turbulence [13]. A stable solid exhibiting het- 
erophase turbulence is a turbulent crystal. The averaging over heterophase configurations, 
described is Section 3, leads to an effective quasiequilibrium picture of heterophase turbulent 
matter, similarly to the averaged treatment of usual turbulence [43] . 

Heterophase turbulence may develop in any condensed matter. Moreover, it is heterophase 
turbulence which is the dynamical cause of all phase transitions. The melting-crystallization 
transition, discussed in Section 4 is not an exception. Heterophase turbulence in a crystal 
shows up as chaotically appearing and disappearing germs of liquid-like regions of disordered 
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lattice. The crystalline and fluid phases can be distinguished between each other by their 
differing symmetries resulting in different features of their spatial local densities (52) and 
(53). Competing phases quite often differ in their symmetry, which, though, is not a manda- 
tory prerequisite for their distinction. Thermodynamic phases may, for instance, having 
the same symmetry, possess different mean densities, as liquid and gas. The possibility of 
heterophase turbulence in a mixture of two coexisting phases with differing densities pi ^ p2 
can be illustrated by a lattice-gas model [44,45]. 

Let the phases be distinguished by their average densities pi and p2, so that 

Pi > P2 . (71) 

Suppose that the system consists of N particles each of which can be located in one of the 
Ng sites of a crystalline lattice {a}, with i — 1, 2, . . . , A^^ and N < Ng. The field operator of 
a particle in a u-phase can be presented as an expansion 

i'A^) = ^inu (^iv <^n(r - aj) (72) 

in 

over the Wannier functions 97n(r — a^), labelled by a multi-index n, with tiy = 0, 1 being the 
occupation variable. Each site can be occupied by not more than one particle, which implies 
the validity of the unipolarity condition 



^inv ^inv 1 • (73) 

n 



The particles are assumed to be well locahzed in the lattice sites, so that the matrix elements 
over the Wannier functions (/7„(r — a^) for the kinetic operator K = — ?i^V^/2mo are 

< mi\k\ni > = 6mn hjKo (74) 

and for the interaction potential $(r — r') are 

< mi,nj\^m'i' ,n'j' > = 5mn' ^if ^nm' Sji' . (75) 

Then the replica Hamiltonian (46) reduces to 

Hy = Wu ^{Kq - p)eiy + ^ij^^iuGju . (76) 

i i^j 

The average phase density (50) takes the form 

i 

The chemical potential p is defined by the normalization condition (51). 
Introducing a quasispin operator 

S- - 1 



(78) 
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transforms the replica Hamiltonian (76) to the Ising-type representation 

H^^E^ + ^wlY. ^ijSiuSj. -w^Yl , (79) 



in which 



The phase density (77) becomes 



^ = ]^E*ii' B, = ix-K^-)^w,^ . (80) 



N 

P. = ^(l + ^.), (81) 



where 

m„ = 



^Y.<S^u> (82) 
* i 

is an analog of a reduced magnetization. Condition (71), owing to Eq. (81), yields 

mi > m2 (pi > P2) • (83) 
Normalization (51) can be rewritten as the equation 

1 1 

n= — = -Wi(l + mi) + -^/;2(l + m2) (84) 

defining the chemical potential = //(n, T) as a function of the parameter (84) and temper- 
ature T. 

The considered two-density lattice model can correspond [13,45] to a stable system whose 
thermodynamic potentials (62) or (63) are lower than the potentials (65) of pure phases, so 
that 

0(r, V, ix) < nl{T, V, fx) , F{T, V, N) < Fl{T, V, N) . (85) 

This means that heterophase turbulence may stabilize condensed matter. That is, turbulent 
crystals should be rather a rule than an exotic exception. Any real solid can exhibit, under 
some conditions, especially in the vicinity of transition points, heterophase turbulence, thus, 
becoming a turbulent crystal. 



6 Ferromagnet with Paramagnetic Fluctuations 

Heterophase fluctuations may arise in systems of different nature. Specific features of the re- 
sulting heterophase matter can be well illustrated by a model of ferromagnet with mesoscopic 
paramagnetic fluctuations [1]. 
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The general form of the rephca Hamiltonian of a lattice spin system reads 

= w,NK, + i wlNU, -wlY, Jij.Si, ■ S,v , (86) 

where is the number of spins; K^, kinetic parameter; U,y, crystalline field; Jj^v, exchange 
interactions; Si^, are spin operators associated with a i^-phase. Ferromagnetic and paramag- 
netic phases can be distinguished by the order parameter 

r]^ = \<Si^>\, (87) 

which, in the case of an ideal periodic lattice, does not depend on the lattice-site index. 
Assume that u = 1 labels the ferromagnetic phase, while i/ = 2, the paramagnetic phase. 
By this assumption, 

m 7^ 1 , m = o. (88) 

The zero order parameter of the paramagnetic phase supposes that the case without external 
magnetic fields is considered. In the presence of the latter, the average magnetization of a 
paramagnet would be nonzero. Then, we could distinguish the phases by the condition 
Vi > V2, similarly to the relations (71) or (83). 

Keeping in mind long-range exchange interactions Jj^v, we shall employ the mean-field 
approximation. The notation 

77 = 771 , w^wi , J =j;^ ^ Jiji (89) 
will be used. Then for the probability of ferromagnetic phase, we find 

which follows from Eq. (43) with 

1 



For the average spin (87), we have 



2?7 = tanh fw^ ^ 77) . (91) 



Here and in the following formulas, we shall set /c^ = 1- It is reasonable to assume that the 

arising magnetic order does not influence much the crystalline properties of matter, so that 
the kinetic terms ATj, and the crystalline fields Un do not depend on the phase index i/, 

Ki^K2, Ui^U2 = A. 

Then, with the notation 

u^j, (92) 
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the ferromagnetic weight (90) reduces to 



= , ■ (93) 

Thus, instead of one equation for the average spin, as it would be in the case of a pure 
phase, we have now two coupled equations, one for the average spin (91) and another for the 

probability of ferromagnetic phase (93). Note that the latter can be treated as an additional 
order parameter that becomes necessary for correctly describing heterophase matter. 

The analysis of the model of the heterophase ferromagnet, in the mean-field approxima- 
tion, shows that there exist two critical temperatures 

The overall behaviour of the system essentially depends on the value of the dimensionless 
parameter (92), which measures the relative intensity of disordering interaction A to ordering 
interaction J. According to the stabihty conditions (44) and (45), paramagnetic fluctuations 
in a ferromagnet may exist only when m > 77^ > 0. Therefore, for some u there can exist the 
nudeation temperature T„, defined by the condition 

w = l (T < T„) , 
w<l {T>T„). 

Depending on the parameter u, we can observe five types of behaviour: 

(i) u<0. 

Pure ferromagnetic phase, with w = 1, is absolutely stable up to the critical temper- 
ature T^, where the system becomes paramagnetic. Heterophase ferromagnet can be only 
metastable, having the critical temperature Tc. 

(ii) 0<u<l 

In the region of temperatures < T < T„, the system is a pure ferromagnet, after which 
paramagnetic fluctuations appear, starting from the nucleation point T„, for T > T„. The 
phase transition ferromagnet-paramagnet becomes a first-order transition, occurring at the 
point To, for which < Tq < T'^. 

(iii) < -u < |. 

The nucleation point T„ exists only for this interval of m, such that < T„ < T^. When 
^ 0, then Tn — > T'^ and if — > 0.5, then — > 0. For > 0.5, there is no nucleation point. 
The first-order phase transition takes place at To- 



(iv) u 



Ferromagnet is always heterophase. The transition temperature Tc becomes a tricritical 
point, where critical indices change by a jump. 

(v) u>l. 

Heterophase ferromagnet is stable. The second-order transition happens at the critical 
temperature Tc. 
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The behaviour of the order parameter, ferromagnetic phase probabihty, and other ther- 
modynamic characteristics is illustrated in the following figures. The order parameter rji, 
defined as the average spin (87), for different values of u, is shown in Fig. 4 as a function of 
0/J, with e = keT. 

The probabihty of ferromagnetic phase Wi = w, satisfying Eq. (93), is presented in Fig. 
5 as a function of the reduced temperature O/J. 
The average magnetization 

M = wri (96) 

as a function of G/J is depicted in Fig. 6. For the values u in the range —0.25 < u < 0, 
there exists a maximum of the magnetization Mm at the temperature T^, where 

Mm = ^ , Tm = "^/^^ . (97) 

4 ' 16arctanh(2V^) 

These quantities are shown in Fig. 7 as functions o{u = A/J. 

The first-order phase transition to the paramagnetic state occurs at temperature, where 
the free energy per particle 

for the heterophasc fcrromagnet, given by the expression 

4cosh (^w^ ^ T]^ , (98) 



f f 2 lA 2 2^1 

w — w-\--]u + wrj — — In 



J V 2J ' J 

coincides with the value 



1 = 1-1'-'- 

corresponding to 77 = and w — 0.5. Functions (98) and (99) are illustrated in Fig. 8 for 
u = 0.6. 

The relative entropy 

^ ~ ^ {E = <H>) (100) 



NT 

takes the form 

1 r / _ 7" \i 

(101) 



-2w^ ^ 77^ + In 



4cosh {^'^ ^ T}^ 



which is depicted in Fig. 9 for different u. 
The specific heat 

c„.ip = r|^ (102) 

in portrayed in Fig. 10, from where it is seen that metastable heterophase ferromagnets, 
with w < 0, possess specific-heat maxima below the critical temperature. 

When the critical point Tc becomes tricritical at u = 3/2, the critical indices change by 
a jump. Thus, for the critical indices a and (3, characterizing the asymptotic behaviour of 
the specific heat and magnetization as 

Cij~(-T)-", M~(-T)^ (103) 
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where 



we obtain 



r = nn -0 , 

T 

c 



In any case, a + 2(3 = 1. If one calculates the critical indices with r tending to a small but 
finite value, instead of r — > —0, then one gets the effective indices 

InCj/ ^ InM 

These are pictured in Fig. 11 for |r| = 10~^. 

Finally, Fig. 12 demonstrates the temperature variation of the average spin, defined by 
Eq. (91), for different values of the parameter (92). Only the solutions for stable states are 
presented. 



7 Conclusion 

Mesoscopic heterophase fluctuations may arise in any kind of condensed matter. This is 
illustrated here for crystals with liquid-like fluctuations, for a two-density lattice mixture, and 
for ferromagnets with paramagnetic fluctuations. Of course, these fluctuations appear not 
always, but their existence depends on the system parameters as well as on thermodynamic 
variables. The presence of such fluctuations can make the system thermodynamically more 
stable. Their appearance may provoke different pretransitional phenomena and transitional 
anomalies. 

The fluctuations are mesoscopic in size, since they represent the germs of competing 
thermodynamic phases, which requires that the characteristic length of such a germ be much 
larger than the mean interparticle distance. At the same time, this characteristic length can 
be much smaller than the size of the whole sample. If the fluctuations are not static, but 
evolve in time, appearing and disappearing, their typical lifetimes are to be essentially longer 
than the local-equilibrium time, but can be much shorter than the observation time, which 
means that such temporal fluctuations are mesoscopic in time. In any case, these fluctuations 
are randomly distributed in space. 

The spatial randomness of the mesoscopic fluctuations makes it possible to realize the 
averaging over heterophase conflgurations and to deflne an effective renormalized Hamilto- 
nian. In the course of this averaging, novel quantities emerge, representing the probabilities 
of thermodynamic phases. These probabilities are deflned as minimizers of the related ther- 
modynamic potentials. 

The theory of heterophase fluctuations [1] provides a general approach for treating a 
wide class of materials whose properties can be strongly influenced by the existence of such 
fluctuations. At the same time, many phenomena occurring in condensed matter, which look 
anomalous without taking account of the mesoscopic fluctuations, can be explained by the 
influence of the latter. 
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Figure Captions 



Fig. 1. The crystalline phase probability as a function of temperature. 
Fig. 2. The liquid phase probability as a function of temperature. 
Fig. 3. The glassy phase probabihty versus temperature. 

Fig. 4. The order parameter rji = rj, given by Eq. (91), for different parameters (92): 
u — 0.5 (curve 1), u — 0.6 (curve 2), u — 0.8 (curve 3), u — 1.5 (curve 4), u — oo (curve 5), 
U — —1 (curve 6),u — —0.25 (curve 7),u — —0.1 (curve 8),u— —0.01 (curve 9),u — —0.001 
(curve 10). 

Fig. 5. The ferromagnetic phase probability wi = w, defined by Eq. (93), as a function 
of relative temperature, for the same values of u as in Fig. 4. 

Fig. 6. The average magnetization (96) as a function of relative temperature for different 
values of u as in Fig. 4. 

Fig. 7. The maximum of magnetization and the related temperature, as defined in 
Eq. (97), versus u = A/J. 

Fig. 8. Free energies per particle, given by Eqs. (98) and (99) for u = 0.6. 

Fig. 9. Relative entropy (101) for the values of u as in Fig. 4. 

Fig. 10. Specific heat (102) for different u: u = 0.6 (curve 2), u = 0.8 (curve 3), u — 1.5 
(curve 4), M = 3 (curve with no number), u = oo (curve 5), u = —1 (curve 6). 

Fig. 11. Effective critical indices (105) for r = — 10~^, versus u = A/ J. 

Fig. 12. Average spin, satisfying Eq. (91), for different values of which are denoted 
by the numbers close to the related curves. 
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